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Abstract 

The Groverian measures are analytically computed in various types of three-qubit states. The 
final results are also expressed in terms of local-unitary invariant quantities in each type. This fact 
reflects the manifest local-unitary invariance of the Groverian measure. It is also shown that the 
analytical expressions for various types have correct limits to other types. For some types (type 4 
and type 5) we failed to compute the analytical expression of the Groverian measure in this paper. 
However, from the consideration of local-unitary invariants we have shown that the Groverian 
measure in type 4 should be independent of the phase factor ip, which appear in the three-qubit 
state This fact with geometric interpretation on the Groverian measure may enable us to 
derive the analytical expressions for general arbitrary three-qubit states in near future. 



PACS numbers: 03.67.Mn, 03.65.Ud, 03.67.Bg 



I. INTRODUCTION 



Recently, much attention is paid to quantum entanglement It is believed in quantum 
information community that entanglement is the physical resource which makes quantum 
computer outperforms classical one 2| . Thus in order to exploit fully this physical resource 
for constructing and developing quantum algorithms it is important to quantify the entan- 
glement. The quantity for the quantification is usually called entanglement measure. 

About decade ago the axioms which entanglement measures should satisfy were studied^. 
The most important property for measure is monotonicity under local operation and classical 



communication(LOCC)[4j]. Following the axioms, many entanglement measures were con 
structed such as relativeVntrowfl.eutanglem^tofdistillationBandformatioJflB" 



geometric measure 
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12 



m 



14| , Schmidt measure 15| and Groverian measure 16| . Entan 



glement measures are used in various branches of quantum mechanics. Especially, recently, 
they are used to try to understand Zamolodchikov's c-theorem[l3] more profoundly. It may 
be an important application of the quantum information techniques to understand the effect 
of renormalization group in field theories [rs|. 

The purpose of this paper is to compute the Groverian measure for various three-qubit 
quantum states. The Groverian measure G(if)) for three-qubit state is defined by G(ip) = 
VI - Pmax where 

P max = max Kgilfelfel^)! 2 - (1-1) 
|gi},|g2>,|?3) 

Thus Pmax can be interpreted as a maximal overlap between the given state and product 



states. Groverian measure is an operational treatment of a geometric measure. Thus, 
if one can compute G(ip), one can also compute the geometric measure of pure state by 
G 2 (ip). Sometimes it is more convenient to re-express Eq. fll.ip in terms of the density 
matrix p = \if))(if)\. This can be easily accomplished by an expression 



max Tr IpR 1 



R 2 <g> R 3 ] 



;i-2) 



where R 1 = \qi)(qi\ density matrix for the product state. Eq. fll.ip and Eq. fll.2l) manifestly 
show that Pmax an d G(ip) are local-unitary(LU) invariant quan tities. Since it is well-known 

,5), 
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2jJ], say Ji(i 



that three-qubit system has five independent LU-invariants 
we would like to focus on the relation of the Groverian measures to LU-invariants J,'s in 
this paper. 
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This paper is organized as follows. In section II we review simple case, i.e. two-qubit 
system. Using Bloch form of the density matrix it is shown in this section that two-qubit 
system has only one independent LU-invariant quantity, say J. It is also shown that Grove- 
rian measure and P ma x for arbitrary two-qubit states can be expressed solely in terms of J. 
In section III we have discussed how to derive LU-invariants in higher-qubit systems. In 
fact, we have derived many LU-invariant quantities using Bloch form of the density matrix 
in three-qubit system. It is shown that all LU-invariants derived can be expressed in terms 



of Jj's discussed in Ref.[20|. Recently, it was shown in Ref. [22] that P ma x for ^-qubit state 



can be computed from (n — l)-qubit reduced mixed state. This theorem was used in Ref. [23] 



and Ref. [24j to compute analytically the geometric measures for various three-qubit states. 
In this section we have discussed the physical reason why this theorem is possible from the 
aspect of LU-invariance. In section IV we have computed the Groverian measures for various 
types of the three-qubit system. The five types we discussed in this section were originally 



developed in Ref. [20] for the classification of the three-qubit states. It has been shown that 
the Groverian measures for type 1, type 2, and type 3 can be analytically computed. We 
have expressed all analytical results in terms of LU-invariants Jj's. For type 4 and type 5 the 
analytical computation seems to be highly nontrivial and may need separate publications. 
Thus the analytical calculation for these types is not presented in this paper. The results of 
this section are summarized in Table I. In section V we have discussed the modified W-like 
state, which has three-independent real parameters. In fact, this state cannot be categorized 
in the five types discussed in section IV. The analytic expressions of the Groverian measure 



for this state was computed recently in Ref. 24j. It was shown that the measure has three 



different expressions depending on the domains of the parameter space. It turned out that 
each expression has its own geometrical meaning. In this section we have re-expressed all ex- 
pressions of the Groverian measure in terms of LU-invariants. In section VI brief conclusion 
is given. 

II. TWO QUBIT: SIMPLE CASE 

In this section we consider P ma x for the two-qubit system. The Groverian measure for two- 
qubit system is already well-known [25]. However, we revisit this issue here to explore how the 



measure is expressed in terms of the LU-invariant quantities. The Schmidt decomposition 26] 
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makes the most general expression of the two-qubit state vector to be simple form 



A |00) + A 1 |11) 



(2.1) 



with Ao, Ai > and Aq + Af = 1. The density matrix for \ip) can be expressed in the Bloch 
form as following: 



p = = - [1 <E> 1 + v la a a ®t+ v 2a t <g> cr a + g a pa a ® op\ , (2.2) 



where 



/ 



Vi = v 2 








( 2A Ai 0^ 



9aP 



V 



-2A Ai 
1 



(2.3) 



In order to discuss the LU transformation we consider first the quantity Ua a U^ where U 
is 2 x 2 unitary matrix. With direct calculation one can prove easily 



(2.4) 



where the explicit expression of O a p is given in appendix A. Since O a p is a real matrix 
satisfying OO t = O t O = 1, it is an element of the rotation group 0(3). Therefore, 
Eq. (12.41) implies that the LU-invariants in the density matrix (12.21) are | -£T 2 1 , Tr^g 1 ] etc. 

All LU-invariant quantities can be written in terms of one quantity, say J = ApAf. In 
fact, J can be expressed in terms of two-qubit concurrence [s| C by C 2 /4. Then it is easy to 
show 



\vi\ 2 = \v 2 \ 2 = 1 -47, 

9ctP9afi = 1 + 8 J. 



(2.5) 



It is well-known that P ma x is simply square of larger Schmidt number in two-qubit case 



P max = max (Ao, A?) . 
It can be re-expressed in terms of reduced density operators 



(2.6) 



1 



1 + v 7 ! - 4det/^ 



(2.7) 



4 



where p A = Tr#p = (1 + vi a a a )/2. Since P ma x is invariant under LU-transformation, it 
should be expressed in terms of LU-invariant quantities. In fact, P ma x in Eg. (12.71) can be 
re-written as 

1 (2-8) 



1 + V1-4J 



Eg. (12.81) implies that P ma x is manifestly LU-invariant. 



III. LOCAL UNITARY INVARIANTS 



The Bloch representation of the 3-qubit density matrix can be written in the form 
1 



P 



1 ® 1 ® 1 + v la a a ® 1L ® 1L + v 2a t ® a a ® 1 + w 3q 1 ® 1 ® a 



(3.1) 



+ ®> da ®> CT/3 + <8> 1 <8> CX/3 + ti£p(Ta <8> CT/3 <E> 1 + Qa^^a ® &/3 <%> 0" 7 



,(2) 



(3), 



where <r a is Pauli matrix. According to Eq. (I2.4[) and appendix A it is easy to show that 
the LU-invariants in the density matrix (13.11) are |"U 2 |, |"U 3 |, Tr[h^h^ 1 ' T ), Tr[h^h^ 2 ' T ], 

Tj[/ l (3)/ l (3)T] ) g apj g a ^ etc. 

Few years ago Acm et al|20j represented the three-qubit arbitrary states in a simple form 
using a generalized Schmidt decomposition 26] as following: 



|V) = A |000) + A^lOO) + A 2 |101) + A 3 |110) + A 4 |lll) (3.2) 
with Aj > 0, < if < 7r, and X^A 2 = 1. The five algebraically independent polynomial 



LU-invariants were also constructed in Ref. 



20j: 



Ji = A X A 4 + A 2 A 3 — 2A1A2A3A4 cos <p, 

J2 = AqA^, J3 = AgAg, J4 - 

J5 — Aq(Ji + A2A3 — A^A 2 ). 



(3-3) 



\ 2 \ 2 



In order to determine how many states have the same values of the invariants J±, J2, ... Js, 
and therefore how many further discrete-valued invariants are needed to specify uniquely a 
pure state of three qubits up to local transformations, one would need to find the number 
of different sets of parameters ip and Aj(z = 0,1, ...4), yielding the same invariants. Once 
Ao is found, other parameters are determined uniquely and therefore we derive an equation 
defining Aq in terms of polynomial invariants. 



(Ji + J 4 )A* - (J 5 + </ 4 )A 2 + J 2 J 3 + J 2 J 4 + J 3 J 4 + J\ = 0. (3.4) 

This equation has at most two positive roots and consequently an additional discrete- 
valued invariant is required to specify uniquely a pure three qubit state. Generally 18 LU- 
invariants, nine of which may be taken to have only discrete values, are needed to determine 
a mixed 2-qubit state [271 ]. 

If one represents the density matrix as a Bloch form like Eg. (13. II) . it is possible to 

construct v± a , v 2a , v 3a , h^l, h^, h^l, and g a ^ explicitly, which are summarized in appendix 
B. Using these explicit expressions one can show directly that all polynomial LU-invariant 
quantities of pure states are expressed in terms of Jj as following: 



N 2 = 1 - 4( J 2 + h + Ji) 
12 -1-A(J 1 + J 2 + J 4 ) 



v 3\ 



|w 2 | 2 = 1 -4(Ji + J3 + J4) 

1 + 4(2 J! - J 2 
l-4(Ja 



(3.5) 



Tl 



J*) 

J2 - 2J 3 ) 



TT[h^h^ T ] 

Tr[h^h^ T ] = 1 - 4(Ji - 2J 2 + J 3 ), Tr[hWh® 
1 + 4(2 J x + 2 J 2 + 2 J 3 + 3 J 4 ) 
l-4(Jx + J 2 + J 3 + J 4 - J 5 ). 

Recently, Ref. [221] has shown that P max for n-qubit pure state can be computed from 



a af3 U a u f3 



(n — l)-qubit reduced mixed state. This is followed from a fact 

R n ] = max Tr \pR l ® R 2 ® • ■ ■ ® R n ^ ® ll (3.6) 



max Tr [pi? 1 ® i? 2 



which is Theorem I of Ref. 22j. Here, we would like to discuss the physical meaning of 
Eq. (13.61) from the aspect of LU-invariance. Eq. (13. 6ft in 3-qubit system reduces to 



max Tr Ip^R 1 ® R 2 ] 



fi'R 2 



where p^ 45 = Tr c p. From Eq. (l3.ip p AB simply reduces to 

1 



P= 4 

,0) 



(3) 



a/3 1 - 



(3.7) 



(3.8) 



where t>i Q , v 2a and faig are explicitly given in appendix B. Of course, the LU-invariant 
quantities of p AB are \vi\, \v 2 \, Tr[h^h^ T ], h^v\ a v 2 p etc, all of which, of course, can be 
re-expressed in terms of J\, J 2 , J 3 , J4 and J5. It is worthwhile noting that we need all Jj's 
to express the LU-invariant quantities of p AB . This means that the reduced state p AB does 
have full information on the LU-invariance of the original pure state p. 



Indeed, any reduced state resulting from a partial trace over a single qubit uniquely 
determines any entanglement measure of original system, given that the initial state is 
pure. Consider an (n — l)-qubit reduced density matrix that can be purified by a single 
qubit reference system. Let be any joint pure state. All other purifications can be 
obtained from the state by LU-transformations U ® l®( n ~ 1 ) ) where U is a local unitary 
matrix acting on single qubit. Since any entanglement measure must be invariant under 
LU-transformations, it must be same for all purifications independently of U. Hence the 
reduced density matrix determines any entanglement measure on the initial pure state. That 
is why we can compute P max of n-qubit pure state from the (n — l)-qubit reduced mixed 
state. 

Generally, the information on the LU-invariance of the original n-qubit state is partly lost 
if we take partial trace twice. In order to show this explicitly let us consider p A = TybP AB 
and p B = TyaP AB ' 

P A = g I 1 + VlaVa) (3.9) 
P B = + V 2a<7a] ■ 

Eq. (12.41) and appendix A imply that their LU-invariant quantities are only \v\\ and l^l 
respectively. Thus, we do not need J5 to express the LU-invariant quantities of p A and 
p B . This fact indicates that the mixed states p A and p B partly loose the information of the 
LU-invariance of the original pure state p. This is why (n — 2)-qubit reduced state cannot 
be used to compute P max of n-qubit pure state. 



IV. CALCULATION OF P r 
A. General Feature 



If we insert the Bloch representation 

B} = 1 + gV " B? = t + g >- & (4.1) 
2 2 v ; 

with = [ 1 = 1 into Eq.f l3.7p . P max for 3-qubit state becomes 

Pmax = 7 max [1 + fi • si + f 2 • s 2 + gijSiiS 2j \ (4.2) 
4 |ai|=|3b|=i 
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where 

f x = Tr [p A a] (4.3) 

f 2 = Tr [p B a] 

g iS = Tr [p AB ai <g> a,] . 

Since in Eq. fl4.2p P max is maximization with constraint = | s* 2 1 = 1, we should use the 
Lagrange multiplier method, which yields a pair of equations 

f 1 + gs 2 = A 1 s 1 (4.4) 
^2 + g T si = A 2 s 2 , 

where the symbol g represents the matrix g^- in Eq. fl4.3p . Thus we should solve s\, s 2 , Ai 
and A 2 by eq. (14.41) and the constraint \si\ = \s 2 \ = 1. Although it is highly nontrivial to 
solve Eq. fl4.4p . sometimes it is not difficult if the given 3-qubit state \ip) has rich symmetries. 
Now, we would like to compute P max for various types of 3-qubit system. 

B. Type 1 (Product States): J x = J 2 = J 3 = J 4 = J 5 = 

In order for all Jj's to be zero we have two cases A = J\ = or A 2 = A 3 = A 4 = 0. 
1. A = J\ = 

If A = 0, in Eq.(JOP becomes \if>) = |1) ® |5C) where 

l-BC) = A ie ^|00) + A 2 |01) + A 3 |10) + A 4 |ll). (4.5) 

Thus P m ax for \ip) equals to that for \BC). Since \BC) is two-qubit state, one can easily 
compute P m ax using Eq.f l2.7l) . which is 



P = I 

1 max 2 



1 



/ l-4det(Tr B |5C)(5C|) =- 1 + . (4.6) 



If, therefore, A = J\ = 0, we have P ma x = 1, which gives a vanishing Groverian measure. 



S 



2. A 2 = A 3 = A 4 = 



In this case \ip) in Eg. (13.21) becomes 



|V>=(A o |0>+Aie^|l>)®|0)®|0). 



(4.7) 



Since \ip) is completely product state, P ma x becomes one. 

C. Type2a (biseparable states) 

In this type we have following three cases. 

1. J\ ^ and J 2 = J3 = J4 = J5 = 

In this case we have Ao = 0. Thus P ma x for this case is exactly same with Eq. (I4.6j) . 

2. J 2 / and J\ = J 3 = J 4 = J 5 = 

In this case we have A 2 = A4 = 0. Thus P ma x for equals to that for \AC), where 



|A<7> = A |00> + Aie^llO) + A 2 |ll). 



(4.8) 



Using Eq. fl2.7l) . therefore, one can easily compute P ma x, which is 



Pmax = ~ 1 + V 1 - 4J 2 



(4.9) 



3. J 3 ^ and J\ = J 2 = J4 = J5 = 



In this case P ma x for \ip) equals to that for \AB), where 



\AB) = Ao|00) + A ie ^|10) + A 3 |ll). 



(4.10) 



Thus P max for |V>) is 




(4.11) 



9 



D. Type2b (generalized GHZ states): J 4 ^ 0, J 4 = J 2 = J 3 = J 5 = 



In this case we have Ai = A2 = A3 = and becomes 



|^) = Ao|000) + A 4 |lll) 



(4.12) 



with An + A4 = 1. Then it is easy to show 



n = Tr[p A a] =(0,0,A^-A2 

f 2 = Tr [p B a] = (0, 0, A 2 , - A, 2 






(4.13) 



/ n n n \ 



gij = Tr [p AB a % ® a 3 ] 



1 



Thus P m nr reduces to 



- max [l + (Aq - Xl){si z + s 2z ) +s lz S2z\ 

4 |si| = |s 2 |=l 



(4.14) 



Since Eq. fl4.14p is simple, we do not need to solve Eq. (I4.4|) for the maximization. If Ao > A4, 
the maximization can be achieved by simply choosing s 1 = 5*2 = (0,0,1). If Ao < A4, we 
choose s 1 = s*2 = (0, 0,-1). Thus we have 



Pmax = max(A , A 4 ). 



(4.15) 



In order to express P ma x in Eq. (14.150 in terms of LU- invariants we follow the following 
procedure. First we note 

1 



max — 2 [(Ao + A4) + I Aq - A4I] 



(4.16) 



Since \\\-\\\ = y/(\% + XI) 2 - 4A§A| = VI - 4J 4 , we get finally 

1 



1 + Jl - 4J 4 



(4.17) 



E. Type3a (tri-Bell states) 



In this case we have Ai = A 4 = and lib) becomes 



Ao|000) + A 2 |101) + A 3 |110) 



(4.18) 
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with A§ + XI + A| = I. If we take LU 
into which is usual W-type state 



transformation a x in the first-qubit, 
2811 as follows: 



= Ao|100) + A 3 |010) + A 2 |001). 
The LU-invariants in this type are 



is changed 



(4.19) 



2 \2 
2 A 3 



J\ — A 2 A 
J 3 



\ 2 \ 2 

A A 3 



■h 



\ 2 \ 2 
A A 2 

2A0A2A3. 



(4.20) 



Then it is easy to derive a relation 



hh + Jih + J2J3 = VJ1J2J3 = \h- (4.21) 



Recently, P m ax for \ip') is computed analytically in Ref. [23] by solving the Lagrange 
multiplier equations (14.41) explicitly. In order to express P max explicitly we first define 

(4.22) 





= ^l + xl 


X 2 




— Ao + Xl 


-Al 


r 3 


= A^ + A^ 


-xl 




— 2A0A3. 





Also we define 

a = max(Ao, A2, A3) (4.23) 
b = mid(A , A 2 , A 3 ) 
c = min(A , A 2 ,A 3 ). 

Then P max is expressed differently in two different regions as follows. If a 2 > b 2 + c 2 , P max 



becomes 



a = max(A , A 2 , A 3 ). 



(4.24) 



In order to express P ma x in terms of LU-invariants we express Eq. fl4.24p differently as 

P> ax = \ [{Xl + Xj + Xl) + \X 2 + Xl- Xt\ + |A§ - A^ + A^| + [AS — A3 — A||] . (4.25) 
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Using equalities 



I An + Aq — \i 



I An Aq + An I 



IX 2 X 2 X 2 I 
l A ~~ A 3 ~~ A 2\ 



1 - 4A 2 A 2 - 4A 2 A 2 



1 - 4A 2 A 2 - 4A 2 A 2 



1 - 4A 2 A 2 



4A 2 A§ 



1 - 4(J X + J 2 ) 
1 - 4( Ji + J 3 ) 
1 - 4(J 2 + J 3 ] 



(4.26) 



we can express P, 
If a 2 < 6 2 + c 2 , P 



„ ax in Eq. (14.241) as follows: 
1 + y/1 - 4( Ji + J 2 ) + v 7 ! - 4( Ji + J 3 
becomes 



1 - 4(J 2 + J 3 



^ ^ cua/ (c<j 2 + r\ — rl)(uj 2 + r\ — r|) — rir 2 r 3 



(4.27) 



(4.28) 



It was shown in Ref. [23j| that P ma x = 4R 2 , where R is a circumradius of the triangle Ao, A 2 
and A 3 . When a 2 < b 2 +c 2 , one can show easily r\ = y/l — 4(J 2 + J3), r 2 = a/1 — 4(Ji + J3), 



r 3 



a/1 — 4(J X + J 2 ), and u; = 2v^3- Using u 2 



8AqA 2 A 2 , One can show 



easily that P maa ; in Eq. (l4.28p in terms of LU-invariants becomes 



P' 



(4.29) 

4( Ji + J 2 + J 3 ) - 1 1 ; 

Let us consider Aq = limit in this type. Then we have J 2 = J 3 = 0. Thus P^ ax reduces 



to (1/2) (1 + \J\ — 4 Ji) which exactly coincides with Eq. (l4.6j) . By same way one can prove 
that Eq. ( 14.271) has correct limits to various other types. 



F. Type3b (extended GHZ states) 



This type consists of 3 types, i.e. \\ = A 2 = 0, Ai = A 3 = and A 2 = A 3 = 0. 



1. Ai = A 2 = 



In this case the state (I3.2p becomes 

|V) = A o |000> + A 3 |110> + A 4 |1H) 
with Aq + A 2 + A 2 = 1. The non-vanishing LU-invariants are 



J?, — A n A 



2 \2 



0^3' 



Ja — AnA/. 



(4.30) 



(4.31) 
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Note that J3 + J4 is expressed in terms of solely Aq as 



J3 + Ja = Ajj(l-Ajj) 



(4.32) 



Eq. fl4.30p can be re- written as 



|V) = A |00g 1 ) + - v /l-A2|llg 2 ) 



(4.33) 

where |gi) = |0) and \q 2 ) = (l/y/l — Ao)(As|0) + A4 1 1 ) ) are normalized one qubit states. 

is 



Thus, from Ref.l23|, P max for 



P max = max (Aq, 1 - Aq) = - 



1 + A /(1-2A^ 



(4.34) 



With an aid of Eq. (14. 32[) P max in Eq. (I4.34p can be easily expressed in terms of LU-invariants 
as following: 



1 

2 L 



1 + sf\ - 4(J 3 + J 4 ) 



(4.35) 



If we take A3 = limit in this type, we have J3 = 0, which makes Eq. (14.351) to be (1/2)(1 + 



\/l — 4J 4 ). This exactly coincides with Eq. (l4.17l) . 



2. Ai = A3 = 



In this case lib) and LU-invariants are 



|V> = A |0 gi 0) + Jl-Ag|lgal> 



(4.36) 



and 



Jl — AnA 



2 \2 



J4 — A0A4 



(4.37) 



where |<?i) = |0), \q 2 ) = (l/y/l - Ag)(A 2 |0) + A 4 |l)), and Ag + A^ + A^ = 1. The same method 
used in the previous subsection easily yields 

1 



1 + VI - 4( J 2 + Ja) 



(4.38) 



One can show that Eq. (14.381) has correct limits to other types. 
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In this case lib) and LU- invariants are 



|V) = ^/1-Al|g 1 00) + A 4 |g 2 ll) (4.39) 

and 

Ji = A1A4, J 4 = A^ (4.40) 

where |g x ) = (l/y/l - A|)(A |0> + Aie^l)), \q 2 ) = |1), and Ag + A? + A4 = 1. It is easy to 
show 

P = i 

- 1 mm 2 



1 + VI - 4(Ji + J 4 )j . (4.41) 
One can show that Eg. (14.411) has correct limits to other types. 



G. Type4a (A 4 = 0) 

In this case the state vector \ip) in Eg. (13.21) reduces to 

|V>) = A |000) + Aie^|100) + A 2 |101) + A 3 |110) (4.42) 

with Aq + Af + A| + A3 = 1. The non- vanishing LU-invariants are 

Ji = A^A^ J 2 = A^A^ (4.43) 
J3 — A0A3 J5 = 2AQA2A3. 

From Eg. (14.431) it is easy to show 

VJ1J2J3 = \h- (4.44) 

The remarkable fact deduced from Eg. ( 14.431) is that the non- vanishing LU-invariants are 
independent of the phase factor ip. This indicates that the Groverian measure for Eg. ()4.42p 
is also independent of 

In order to compute P ma x analytically in this type, we should solve the Lagrange multiplier 
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equations ()4.4j) with 



r*i = Tr[p A a] = (2AoAi cosy?, 2AoAi sin (p, 2Aq — 1) 
r 2 = Tr[p B a] = (2AiA 3 cos<p, — 2AiA 3 sin<p, 1 — 2Ag) 

7 2A A 3 

-2A A 3 



(4.45) 



gij = Tr \p AB Oi <g> a 



2A Ai cos<p 
2A Ai sin cp 

y-2AiA 3 cos<p 2AiA 3 sin<p Ag - A^ - X\ + A 3 J 



Although we have freedom to choose the phase factor (p, it is impossible to find singular values 
of the matrix g, which makes it formidable task to solve Eq. fl4.4p . Based on Ref.j^ and 



Ref. 



241 ] . furthermore, we can conjecture that P m ax for this type may have several different 



expressions depending on the domains in parameter space. Therefore, it may need long 
calculation to compute P max analytically. We would like to leave this issue for our future 
research work and the explicit expressions of P ma x are not presented in this paper. 



H. Type4b 

This type consists of the 2 cases, i.e. A2 = and A 3 = 0. 



1. A 2 = 

In this case the state vector \ip) in Eq. (13.21) reduces to 

|V>) = A |000) + Aie^|100) + A 3 |110) + A 4 |lll) (4.46) 

with Aq + Xf + A| + A4 = 1. The LU- invariants are 

Ji = Xi\l Js = X 2 Xl Ja = XIX\. (4.47) 

Eq. (14.471) implies that the Groverian measure for Eq. (l4.46p is independent of the phase factor 
p like type 4a. This fact may drastically reduce the calculation procedure for solving the 
Lagrange multiplier equation (14.41) . In spite of this fact, however, solving Eq. (14.41) is highly 
non-trivial as we commented in the previous type. The explicit expressions of the Groverian 
measure are not presented in this paper and we hope to present them elsewhere in the near 
future. 
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2. A 3 = 



In this case the state vector \ip) in Eg. (13.21) reduces to 

|V>) = A |000) + Aie^llOO) + A 2 |101) + A 4 |lll) 
with Ag + Af + A| + A| = 1. The LU- invariants are 



Ji — A^A^ J 2 — AqA^ </ 4 — AqA 



2 \2 



>2\2 



(4.48) 



(4.49) 



Eg. (14.491) implies that the Groverian measure for Eg. (l4.48p is independent of the phase factor 
(p like type 4a. 



I. Type4c (Ai = 0) 



A 2 A 3 J 2 



■■> - A0A3 



In this case the state vector \ip) in Eg. (13.21) reduces to 

|V) = A |000) + A 2 |101) + A 3 |H0) + A 4 |lll) 
with Ag + A 2 + A| + A4 = 1. The LU- invariants in this type are 

Ji 

From Eg. (14.511) it is easy to show 

W2 + J3 + Ji) + hh = = \h- 



,2\2 



AqA 2 J; 



A A 4 J5 



2AQA2A3. 



In this type fi, r 2 and gij defined in Eg. (14.31) are 

ri = (0,0,2Ag-l) 

f 2 = (2A 2 A 4 , 0, Aq + A 2 - A3 - A4) 



/ 



9ij 



\ 





-2A0A3 









l-2Xj 



(4.50) 



(4.51) 



(4.52) 



(4.53) 



2A0A3 


— 2A 2 A 4 

Like type 4a and type 4b solving Eg. (l4.4p is highly non-trivial mainly due to non- 
diagonalization of g^. Of course, the fact that the first component of r 2 is non-zero makes 
hard to solve Eg. (14.41) too. The explicit expressions of the Groverian measure in this type 
are not given in this paper. 
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J. Type5 (real states): tp = 0, n 



1. <p = 



In this case the state vector in Eg. (13.21) reduces to 



A |000) + Ai|100) + A 2 |101) + A 3 |110) + A 4 |H1> 



with Aq + A 2 + A 2 . + A3 + \\ = 1. The LU- invariants in this case are 



J\ — (A2A3 — AlA4) 2 J2 — AqA^ J3 — AgAg 

J4 = AqA 2 J5 = 2AgA2A3(A2A3 — A1A4). 



It is easy to show yj J1J2J3 = t/5/2. 



(4.54) 



(4.55) 



2. if = TT 



In this case the state vector \ip) in Eg. (13.21) reduces to 



A |000) - Ai|100) + A 2 |101) + A 3 |110) + A 4 |lll) 



with Aq + \\ + \\ + A3 + A| = 1. The LU- invariants in this case are 



J\ — (A2A3 + A1A4) 2 J2 — AgAg J3 

J4 = A0A4 J5 = 2AgA2A3(A2A3 + A1A4). 



x 2 x 2 

A A 3 



(4.56) 



(4.57) 



It is easy to show y/ J1J2J3 = J5/2 in this type. 

The analytic calculation of P m ax i n type 5 is most difficult problem. In addition, we don't 
know whether it is mathematically possible or not. However, the geometric interpretation 



of Pmax presented in Ref. 



231 ] and Ref.[24| may provide us valuable insight. We hope to leave 



this issue for our future research work too. The results in this section is summarized in 
Table I. 
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Type 


conditions 


P 

TYltlX 


Type I 


7=0 


1 






j- — n pvpp"nf T-> 


1/1 4. A /l _ 4 7^ 


TVnp TT 
iype 11 




= D excent .To 


A ( 1 4 ./I _ 4J 9 ) 






= o except J3 


A (1 4. /I _ 4 J ) 
2 V v y 




b 


Ji = except J 4 


1 (1 + y/\-UA 

2 \ v 




a 


Aj = A 4 = 


A fl+ A /l-4(J 1 + J 2 )+ v /l-4(J 1 4J3)4 v /l-4(J 2 4J3)) if CL 2 > b 2 + (? 








WJ1J2J3/ (4(Ji + J 2 4- J 3 ) - 1) if < 6^ + 


Type III 




A x = A 2 = 


1(1 + 71-4(73 + 74)) 




b 


Ai = A 3 = 


i ( 1 + v/1 - 4( J 2 + J 4 ) ) 






A 2 = A 3 = 


A ( 1 4. /1 _ 4( j 4- J 4 ) ) 




a 


A, — n 


lTinP'nPTlHPTlf" OT C/T Tinf nrpcprifpn 
111UC U C11UC11 L Ul l±Wu UlCoClltCU. 


Type IV 


b 




lTinpriPTinpnt of in m not nrpcipntpn 






At = 

/V3 w 


indpDPTidprit of go* not nrpspntpd 




c 


Ai = 


not presented 


Type V 


^ = 


not presented 






(f = 7T 


not presented 



Table I: Summary of P max in various types. 



V. NEW TYPE 
A. standard form 

In this section we consider new type in 3-qubit states. The type we consider is 

|$) = a|100) + 6|010) + c|001) + g|lll), a 2 + b 2 + c 2 + q 2 = 1. (5.1) 

First, we would like to derive the standard form like Eq. fl3.2p from |<&). This can be achieved 
as following. First, we consider LU-transformation of |$), i.e. (U <S> 1 <S> 1)1$), where 



U = ^—j- \ V * \ (5-2) 

Vaq + bc I _ v /fc faq 
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After LU-transformation, we perform Schmidt decomposition following Ref . [20] . Finally we 
choose 6 to make all \ to be positive. Then we can derive the standard form (13.21) from |$) 
with ip = or 7r, and 



Ao= / (00 + ^(06 + eg) 
aq + be 

^Jabcq I 2 I 2 r 2 21 

Ai = — ^= |a + g —o — c| 

y (afc + eg) (ac + 6g) (ag + 6c) 

A 2 = t-|oc - bq\ 
Ao 

-^3 = T-\ab - cq\ 
^0 



A 



2^/abcq 



A 



Ao 

,2 1 1,2 1 ^2 1 „2 



It is easy to prove that the normalization condition a + b + c + q — 1 guarantees the 
normalization 

+ A? + Al + A| + Al = 1. (5.4) 

Since |<3>) has three free parameters, we need one more constraint between Aj's. This addi- 
tional constraint can be derived by trial and error. The explicit expression for this additional 
relation is 



KM + ^ + AD = \ - |(A^ + A4XA3 + Xi). (5.5) 



Since all Aj's are not vanishing but there are only three free parameters, |$) is not involved 
in the types discussed in the previous section. 



B. LU-invariants 

Using Eq. (15.31) it is easy to derive LU-invariants which are 

J x = (A X A 4 - A 2 A 3 ) 2 = {ab + cq)2 \ ac + bq) 2 ( 5 - 6 ) 

x \2abcq\a 2 + q 2 — b 2 — c 2 | — (aq + bc)\(ab — cq)(ac — bq)\j 2 
J 2 = A 2 A 2 = (ac - bq) 2 
J 3 = Xl\l = (ab-cq) 2 
J4 = AqA^ = 4a6cg 
J5 = Aq (Ji + A2A3 — A 2 A 2 ) . 



19 



One can show directly that J5 = 2-y/Ji J2J3. Since |$) has three free parameters, there 
should exist additional relation between Jj's. However, the explicit expression may be hardly 
derived. In principle, this constraint can be derived as following. First, we express the co- 
efficients a, b, c, and q in terms of J\, J2, J3 and J4 using first four equations of Eq. (15.61) . 
Then the normalization condition a 2 + b 2 + c 2 + q 2 = 1 gives explicit expression of this addi- 
tional constraint. Since, however, this procedure requires the solutions of quartic equation, 
it seems to be hard to derive it explicitly. 

Since J\ contains absolute value, it is dependent on the regions in the parameter space. 
Direct calculation shows that J\ is 



Ji 



(aq — be) 2 when (a 2 + q 2 — b 2 — c 2 )(ab — cq)(ac — bq) > 
(aq — be) 2 [1 + 2(ab — cq)(ac — bq)(aq + be)/ [ah + cq)(ac + bq)(aq — be)} 2 (5.7) 
when (a 2 + q 2 — b 2 — c 2 )(ab — cq)(ac — bq) < 0. 



Since P ma x is manifestly LU-invariant quantity, it is obvious that it also depends on the 
regions on the parameter space. 



C. calculation of P„ 



Pmax for state |<3>) in Eq.f l5.ll) has been analytically computed recently in Ref. [2J]. It turns 
out that Pmax is differently expressed in three distinct ranges of definition in parameter space. 
The final expressions can be interpreted geometrically as discussed in Ref. |24j] . To express 
Pmax explicitly we define 



7 2 1 2 2 2 2,2 1,2 2 

T\ = c > + c ■ — a — q r 2 = a + c — b — q 

r 3 = a 2 + b 2 — c 2 — q 2 to = ab + qc fi = ab — qc. 



(5.8) 



The first expression of P m ax, which can be expressed in terms of circumradius of convex 
quadrangle is 

p(Q) = 4:{ab + gc){ac + qb)(aq + bc) ^ 



Auj 2 — r 2 



The second expression of P ma x > which can be expressed in terms of circumradius of crossed- 
quadrangle is 



p(GQ) 



(ab — cq)(ac — bq)(bc — aq) 



AS 2 



(5.10) 
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where 

Si = h( a + b + c + q)(a + b-c-q)(a-b + c-q)(-a + b + c-q). (5.11) 
Id 

The final expression of P m ax corresponds to the largest coefficient: 

= max(a 2 ,& 2 ,c 2 ,g 2 ) = \(l + N + |r 2 | + |r 3 |) . (5.12) 



The applicable domain for each P max is fully discussed in Ref . [24 ] . 

Now we would like to express all expressions of P max i n terms of LU-invariants. For the 
simplicity we choose a simplified case, that is (a 2 + q 2 — b 2 — c 2 )(ab — cq)(ac — bq) > 0. Then 
it is easy to derive 

r\ = 1 -4(J 2 + J 3 + J 4 ) r 2 = 1 - 4( + J 3 + J 4 ) (5.13) 
r 2 = 1 - 4( Ji + J 2 + J 4 ) cu 2 = J 3 + J 4 . 

Then it is simple to express Pmax and Pmax^ as following: 



p(Q) = V(Ji + J 4 )(J 2 + J 4 )(J 3 + J 4 ) , , 

4(J X + J 2 + J 3 + 2J 4 ) - 1 1 • J 

p(CQ) 4 a/ J1J2J3 



4(Ji + J 2 + J 3 + Ji) ~ 1' 
If we take q = limit, we have A 4 = J 4 = 0. Thus Pmax and Pmax reduce to 



4a/JiJ 2 J 3 /(4(Ji + J 2 + ^3) - 1), which exactly coincides with P< ax in Eq. (l4.29j) . Finally 
Eq. fl5.13l) makes Pmax to be 

^1 = \ (l + a/1 - 4( J 2 + J 3 + J 4 ) + a/1 - 4( Ji + h + J 4 ) + a/1 - 4(Ji + J 2 + J 4 )) • 

(5.15) 

One can show that Pmax equals to P^ ax in Eq. (I4.27| ) when q = 0. This indicates that our 
results ( I5.14p and (15 . 15[) have correct limits to other types of three-qubit system. 



VI. CONCLUSION 

We tried to compute the Groverian measure analytically in the various types of three- 
qubit system. The types we considered in this paper are given in Ref. [3] for the classification 
of the three-qubit system. 
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For type 1, type 2 and type 3 the Groverian measures are analytically computed. All 
results, furthermore, can be represented in terms of LU-invariant quantities. This reflects 
the manifest LU-invariance of the Groverian measure. 

For type 4 and type 5 we could not derive the analytical expressions of the measures 
because the Lagrange multiplier equations (I4.4p is highly difficult to solve. However, the 
consideration of LU-invariants indicates that the Groverian measure in type 4 should be 
independent of the phase factor (p. We expect that this fact may drastically simplify the 
calculational procedure for obtaining the analytical results of the measure in type 4. The 
derivation in type 5 is most difficult problem. However, it might be possible to get valuable 
insight from the geometric interpretation of P m axi presented in Ref.j^ and Ref. 24]. We 
would like to revisit type 4 and type 5 in the near future. 

We think that the most important problem in the research of entanglement is to under- 
stand the general properties of entanglement measures in arbitrary qubit systems. In order 
to explore this issue we would like to extend, as a next step, our calculation to four-qubit 
states. In addition, the Groverian measure for four-qubit pure state is related to that for 
two-qubit mixed state via purification 29] . Although general theory for entanglement is far 
from complete understanding at present stage, we would like to go toward this direction in 
the future. 

Acknowledgement: This work was supported by the Kyungnam University Research 
Fund, 2007. 
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Appendix A 

One can easily show that the elements of O defined in Eg. (I2.4j) are given by 





= ^ (U U U* 22 + U* n U 2 2 + U 12 U* 21 + U* 12 U 2 l) 


022 


= 2 ( M n M 22 + u *n u 22 - u 12 u* 21 


- u* n u 2 x) 


033 


= Pill - K2I 




012 


1 

= 2 ( M 12 M 21 + M 11 M 22 _ M 12 M 21 


- u* u u 22 ) 


021 


1 

= ^ ( U 12U 21 + U* U U 22 - U* 12 U 2 i 


- u u u* 22 ) 


013 


= u n u* l2 + u* n u 12 




031 


= u u u* 21 + u* n u 21 




023 


= -i (wu«*2 + u* 21 u 22 ) 




032 


= i (u u u 21 + u* 12 u 22 ) 





where Uij is element of the unitary matrix defined in Eq. fl2.4p . It is easy to prove OO t 
O t O = 1, which indicates that O a p is an element of 0(3). 
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Appendix B 

If the density matrix associated from the pure state \ip) in Eg. (13.21) is represented by 
Bloch form like Eg. (13. II) . the explicit expressions for are 



Vl 



( 



\ 



2AoAi COS<y? 

2A Ai siny? 

\ 2 \ 2 X 2 \' 
A — A x - A 2 - A. 



A2 



^2 



' 2AiA 3 cosv2 + 2A 2 A 4 ^ 
— 2AiA 3 sin cp 

\ A + A l + A 2 - A 3 ~ A 4 y 



?'3 



2AiA 2 cos + 2A3A4 
— 2A1A2 sin <p 
\ A o + A i - A 2 + A ; 



\ 



A 4 



and the components of are 



(i> 
11 



h (1) 

' 4 33 



13 

(1) 



2A2A3 + 2A1A4 cos <p, 

A + A l — A 2 — A 3 + A 4' 
— 2A2A4 + 2A1A3 COS (p, 



h 



22 = 2A2A3 — 2AiA4Cos<^ 



'12 



h 



(i> 
21 



-2AiA 4 sin ip 



''31 



-2A3A4 + 2A1A2 cos <p 



h 



h 



v 23 

(2) 
11 

(2) 
12 

( 2 ) 
31 

(2) 
32 



— 2A1A3 sin<£>, 



h%2 = — 2AiA2sin<yj 



-h 



22 



2A A 



2- 



( 2 ) 



^3 



A 



A l + A 2 — A 3 + K 



h 



(2) 
21 



0. 



hf 3 = 2AqAi cos p> 



— 2A3A4 — 2AiA 2 cos <p, 
2A1A2 sin tp. 



(2) 

h 23 = 2A Ai sinyj 



(B.l) 



(B.2) 



(3) (2) 

The matrix h a * is obtained from h a l by exchanging A2 with A3. The non- vanishing compo- 
nents of g al 3-y are 



#111 — — #122 — — #212 — ~ #221 — 2A0A4 

#113 = — #223 = 2A0A3, gisi = —#232 = 2A0A2 

#133 = 2A Ai cos <p, g 23 3 = 2A Ai sin (p 

#312 = #321 = 2AiA 4 sin (p, g 311 = -2A 2 A 3 - 2AiA 4 cos (p 

#313 = 2A 2 A 4 - 2Ai A 3 cos ip, g 32 2 = -2A 2 A 3 + 2AiA 4 cos tp 

#323 = 2AiA 3 sin cp, g 331 = 2A 3 A 4 - 2AiA 2 cos ip 

#332 = 2AiA 2 sin tp, g 333 = Aq - A 2 + A2 + A3 - \\. 



(B.3) 
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